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The epitaxial growth of ultrathin Fe film on Si(l 1 1) surface provides an excellent opportunity to investigate 
the contribution of magnetic anisotropy to magnetic behavior. Here, we present the anisotropic 
magnetoresistance (AMR) effect of Fe single crystal film on vicinal Si(lll) substrate with atomically flat 
ultrathin p(2 X 2) iron silicide as buffer layer. Owing to the tiny misorientation from Fe(lll) plane, the 
symmetry of magnetocrystalline anisotropy energy changes from the six- fold to a superposition of six- fold, 
four- fold and a weakly uniaxial contribution. Furthermore, the magnitudes of various magnetic anisotropy 
constants were derived from torque curves on the basis of AMR results. Our work suggests that AMR 
measurements can be employed to figure out precisely the contributions of various magnetic anisotropy 
constants. 

Magnetic anisotropy is not only the origin of long-range-magnetic-order in low dimensional system 1 , but 
also plays a vital role in determining the magnetic properties for magnetically hard, magnetically soft 2 , 
high-frequency magnetic materials 3 , ultrahigh density magnetic recording media and spintronic mate- 
rials. Recently, the growth and magnetic properties of single crystal Fe film on Si(lll) surface have been 
investigated owing to its application in integration of magnetic devices in Si-based technology and new oppor- 
tunities in spintronics 4 " 11 . The research of Fe ultrathin film magnetization has approached down to the atomic 
level by a powerful spin-polarized scanning tunneling microscopy (SP-STM) 12 " 14 . In the case of bcc Fe film grown 
on Si(lll) substrate, the six-fold symmetry of magnetic anisotropy energy exists only when magnetization is 
confined strictly in the Fe(lll) plane. A small structural modification is sufficient to destroy the six-fold 
symmetry as a result of the contributions from other magnetic anisotropy energies 7,8 . In our previous work, 
we observed that the six-fold symmetry of the in-plane resonance field for Fe(lll) film was changed into the 
superposition of a four-fold and a two-fold contribution due to the presence of atomic step of the vicinal 
substrate 15 . Furthermore, we also observed some difference between ferromagnetic resonance (FMR) results 
and magnetization measurements. FMR results demonstrated that the azimuthal angular dependence of in-plane 
resonance field is a six-fold symmetry with a weak uniaxial contribution, while the remanence of hysteresis loops 
displays a two-fold one 16 . Therefore, the analysis of the magnetic anisotropy and magnetization reversal should be 
carried out carefully for Fe(lll) films on Si(lll) substrate. 

Up to date, various methods, such as magnetic hysteresis loop measurement, torque measurement 17 , ferromag- 
netic resonance (FMR) 7 , rotational magneto-optic Kerr effect (ROT-MOKE) 18 , and magnetic transverse biased 
initial inverse susceptibility and torque (TBIIST) 19 , have been developed to determine the magnetic anisotropy 
constants. Since the coherent domain rotation magnetization reversal for ultrathin film is not always occurred, 
especially when the applied field is lower than saturation field, the detailed information regarding the magnetic 
anisotropy cannot be distinguished precisely from the magnetization hysteresis loops. Alternatively, magnetoran- 
sport method has been proved to be an ideal probe of magnetic anisotropy constants in the thin single layer films 
by anisotropic magnetoresistance (AMR) 20 " 26 , which determines the anisotropy field strength by realization of a 
coherent magnetization reversal (Stoner-Wohlfarth-like). This can be achieved by applying a sufficiently large 
field to guarantee a true single-domain rotation. Here, we carried out the AMR measurements in ultrathin single 
crystalline Fe film on vicinal Si(l 1 1) substrate. On the basis of AMR curves, the angle between the magnetization 
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and magnetic field, and hence the normalized magnetic torque can be 
derived. Finally, the uniaxial magnetic anisotropy, first- and second- 
order magnetocrystalline anisotropy constants were precisely 
obtained by fitting the normalized magnetic toque curves. Our work 
suggests that the extremely sensitive AMR can provide the detailed 
contributions of various magnetic anisotropy constants, including 
the first- and second- order magnetocrystalline anisotropy constants, 
as well as step-induced in-plane uniaxial magnetic anisotropy con- 
stant, of ultrathin Fe single crystal film on vicinal Si(lll) surface. 

Results 

Figure 1(a) shows the schematic configuration of the sample and the 
coordinate system used in our AMR measurements and data ana- 
lysis. Although the substrate supplier declares that the Si(lll) sub- 
strate with orientation accuracy is 0.10° (nominal miscut angle /? = 
0.10°), a local miscut varies from point to point on the Si surface will 
take place due to cutting imperfection or mass transport under the 
direct heating. Recently, we adopted a novel method to tune the 
terrace width of Si(lll) substrate by varying the direction of heating 
current 27 . Large scale images (850 nm X 850 nm) of the Si (111) 
substrate were employed to determine the local variation of miscut 
angles. The typical large scale STM images indicates that the nar- 
rower terraces are companied by a very broad (>400 nm) terrace 
Figure 1(b). From the section line profile along the perpendicular 



direction to the terrace steps (Inset of Fig. 1(b)), a single atomic step 
in Si(l 1 1) surface is about 0.30 nm high can be estimated, which is in 
good agreement with the value reported by Lin et al 28 . On the basis of 
the relation between atomic height h and terrace width w, ft = 
arctan(/z/w), the local miscut angle /? can be various from 0.04° to 
0.30° with a mean miscut angle of 0.10°. The sharp LEED pattern 
plotted in figure 1(b) demonstrates an atomically flat Si(lll)-7 X 7 
recostructured surface. 

Figure 1(c) illustrates the STM image of the iron silicide template 
on the Si (1 1 1) substrate and the corresponding 2X2 LEED pattern. 
The iron silicide template comprises of steps separated by the flat p(2 
X 2) reconstructed terraces. Compared with Si(lll)- 7 X 7 recos- 
tructured surface, the step edges in the STM image for p(2 X 2) iron 
silicide reconstructed surface are not so sharp owing to the random 
diffusion of Fe atoms on Si substrate and the intermixing between Fe 
and Si atoms. The atomically flat terraces are generally used as a 
template for preparing ultrathin Fe single crystalline film. Fig. 1(d) 
shows the STM image of the iron deposited on p(2 X 2) iron silicide 
(lll)/Si(lll) surface for 21 ML and the LEED pattern. The LEED 
pattern indicates that three-fold symmetry still exists even for 
a thickness reaching 21 ML, suggesting a bcc Fe(lll) film. The 
epitaxial relationships between the Fe(lll) film, the iron sili- 
cide template and the Si substrate are following 29 : Fe(lll)||p 
(2 X 2)(lll)||Si(lll) and Fe[-1-12] ||p(2 X 2) iron silicide 
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Figure 1 | The schematic configuration of the sample preparation and coordinate system used in AMR measurements and data analysis, and STM 
images of sample surface, (a) Schematic configuration of the sample and the coordinate system used in our AMR measurements and data analysis; 
(b) STM image and LEED patterns of 0.1° vicinal Si(l 1 1)-7 X 7 reconstructured surface; (c) STM image and LEED patterns of 2 X 2 iron silicide phases 
obtained at annealing temperature T a = 700 K for 10 min after 1.5 ML Fe deposited on 0. 1 ° vicinal Si( 1 1 1 ) surface; (d) STM image and LEED patterns of 
21 ML Fe onp(2 X 2)/Si(lll). 
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[-1-12] ||Si[ll-2]. Owing to the large lattice mismatch, STM topo- 
graphic image shows that Fe atoms aggregate immediately into three- 
dimensional (3D) islands. The first stage is the growth of uniformly 
strained wetting layers, and elastic energy increases fast with film 
thickness. The constrained film becomes unstable at a critical thick- 
ness and three dimensional islands appear for strain relief. With 
further increasing film thickness, misfit dislocations are introduced 
to further relieve stain resulting in deeper facets, dome and ridge- 
trough structures. The section line in the insert of figure 1(d) indi- 
cates that the surface corrugation of grainy thin Fe(lll) is rather 
large. In our previous work, the effect of strain at the Fe/FeSi interface 
on the magnetic anisotropy has been discussed 30 . Since thickness of 
Fe film (about 21 ML) is far thicker than the critical thickness, the 
strain is released, and consequently has no significant influence on 
magnetic anisotropy constants. 

The total free energy density of the system with the external field H 
is considered as the following formula 7 : 

E= — [i 0 M s H'ju + Ki ((x\(x\ + (x\tx\ + tx\{x\) +K 2 ((x 2 l (xl(xl) 

-K u {fi-h u f +K d (fi-h d f -K s {fi-h s f 

where the first term is the Zeeman energy, fi is the unit vector of the 
magnetic vector and M s is the saturation magnetization of Fe (taken as 
the bulk value 1.74 X 10 6 A/m); the second and third terms are cubic 
magnetocrystalline anisotropy energy, a,- represents the directional 
cosines of the magnetic vector fi with respect to the cubic axes 
[100], [010] and [001], K Y and K 2 are the first two cubic magneto- 
crystalline anisotropy constants; the last three terms sequentially refer 
to the uniaxial magnetic anisotropy energy, and surface magnetic 
anisotropy energy and out-of-plane demagnetization energy. K w K d 
and K s are the corresponding magnetic anisotropy constants. The unit 
vector h u with its orientation along the step direction represents the 
direction of the easy axis of the uniaxial magnetic anisotropy. h d and 
n s are the unit vectors normal to vicinal (111) film plane and the (111) 
plane, respectively. It should be noted that the unit vector h d is per- 
pendicular to the vicinal plane, which is different from a simple flat 
thin (111) crystal plane with its hard axis of the out-of-plane. 

Figures 2 (a) and (b) present the angular dependence of the first- 
and second-order magnetocrystalline anisotropy energy terms in the 
Fe(lll) plane along [11-2] with various miscut angles, where K Y = 
4.5 X 10 4 J/m 3 and K 2 = O.OSK^ 6 , respectively. We can find that the 
Ki energy term is invariable (solid line circle in Fig. 2(a)) and the K 2 
energy term is six-fold symmetry in exact Fe(lll) plane, i.e. miscut 
angle /? = 0°. However, the K x energy term can be changed to a four- 
fold symmetry by a slight misorientation from (111) plane, i.e. /? 7^ 
0°. Figure 2(b) demonstrates that the symmetry of the K 2 energy term 
keeps unchanged. 

Since the magnetization reversal process is largely governed by the 
symmetry, magnitudes and directions of the competing magnetic 



anisotropy energies, the symmetry of magnetic anisotropy energy is 
usually probed by magnetic hysteresis loop. The MOKE hysteresis 
loops at various angles cp H between the [-110] axis and magnetic field 
H, indicate that the easy axis is perpendicular to the step direction, (p H 
= 90° (Fig. 3(a)). Similar phenomena have been reported in the system 
Fe/W(001) or Au/Co/Cu/Si(lll), which have been explained by the 
step-induced anisotropy 31 ' 32 . Unfortunately, owing to the small coer- 
civity (<10 Oe) in the sample, only the two-fold symmetry in reman- 
ence and coercivity can be confirmed from figures 3(b) and 3(c), 
respectively. The contribution of magnetocrystalline anisotropy con- 
stants Ki and K 2 cannot be determined from MOKE measurements. 

In order to figure out the magnetocrystalline anisotropy constants 
Ki and K 2 , we carried out the AMR measurements. We found that the 
resistances of Si substrate with and without iron silicide buffer layer, 
which are almost the same values, are quite larger than the resistance 
of Fe ultrathin film. Furthermore, the Si substrate and iron silicide 
buffer layer have no contribution to AMR. Therefore, the metallic Fe 
single crystal film gown on iron silicide buffer layer and Si(lll) 
surface provides an ideal system to perform AMR measurements. 
In the case of Fe single- crystalline system, the AMR can be expressed 
as 20 " 26 : 

Rxx =R± + {R\\-R±) cos 2 cp M (2) 

where (p M is the angle between the Fe magnetization M Fe and the 
current flow I, R\\ and R± is the resistance at cp M = 0° and cp M = 90°, 
respectively. 

Figure 4(a) shows the angular dependence of the in-plane AMR 
with different applied fields. The external magnetic fields are larger 
than saturation field to guarantee a true single-domain rotation and 
eliminate the ordinary magnetoresistance effect. During rotation of 
the sample, the AMR values show an oscillated behavior between the 
maximum value R\\ and minimum value R± 9 respectively. However, 
owing to the magnetic anisotropy, M Fe is no longer kept along with the 
external field H during rotation, i.e. (p M < (p H . Therefore, the AMR 
curves do not follow the cos 2 (p H relationship. The correlation between 
(p H and (p M can be obtained from Fig. 4 (a) and plotted in Fig. 4 (b). 

On the basis of the angle difference between (p H and (p M , we can 
further calculate the magnetic torque L((p M ) = /n 0 M s H sin (cp H — <p M ) 
curves from Fig. 4(b) at different external fields. In order to compare 
magnetic torques at different fields, the normalized magnetic torque 
K<Pm) = L(q> M ) I ^ 0 M S H = sin (g? H — cp M ) was introduced. As shown 
in figure 5(a), the normalized magnetic torque curves exhibit differ- 
ent shapes with different external field H. In equilibrium state, the 
torque acting on M Fe due to His equal in magnitude to the torque due 
to the magnetic anisotropics of the sample. Since the demagnetiza- 
tion field is normal to the Si(lll) plane, its contribution to the 
magnetic torque is zero. According to Eq.(l), the normalized mag- 
netic torque can be written as: 
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Figure 3 | The typical hysteresis loops and angular dependence of in- 
plane magnetic coercivity and remanence. Experimental hysteresis loops 
at 0 and 90 deg. (a), angular dependence of in-plane magnetic coercivity 
(b) and remanence (c) by means of MOKE. 
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Figure 4 | The in-plane AMR curves at different curves and the 
correlation between p H and Pm- Angular dependence of the in-plane AMR 
(a), and the correlation between (p H and (p M at different fields (b). 

where y = -arccos ^1 j V?^J + /?,/? is the miscut angle of the substrate. 

Although the value of iC s (~10 6 J/m 3 ) is far larger than that of 
iC u (~10 2 J/m 3 ) for ultrathin Fe film 1516 , we can calculate from 
Eq.(3) that the value of torque contributed by K s is at least two order 
of magnitude smaller than that contributed by K u . Therefore, the 
contribution of surface anisotropy constant to the torque can be 
neglected. 

It can obviously from Eq.(4) that the magnetic torque shows a 
superposition of two-, four- and six-fold magnetic anisotropics from 
the step -induced uniaxial magnetic anisotropy K u , the first- order 
magnetocrystalline constant K Y and the second-order magnetocrys- 
talline anisotropy constant K 2 , respectively. The two-fold symmetry 
disappears gradually with increasing external field H, suggesting that 
the strength of K u is very weak. Therefore, in order to distinguish the 
contribution from K u , the external field H should be kept slightly 
larger than the saturation field. 

From Eq.(3), we can obviously find that the normalized torque 
l((p M ) is significantly affected by the substrate's miscut angle The 
tendency of the anisotropy energy is complicated. We can find that 
the four-fold anisotropy energy changes significantly (Fig. 2(a)), 
while the six-fold anisotropy energy almost does not change with 
the miscut angle (Fig. 2(b)). In the case of /? = 0.0°, the K x term is 
zero. Usually the miscut angle of the substrate cannot be neglected, 
and thus the contribution from the first- order magnetocrystalline 
anisotropy constant in vicinal (111) plane must be taken into 
account. 

In order to investigate the tiny variation of miscut angles on the 
fitting parameters, Figure 6 illustrates the fitted magnetic anisotropy 
constants for various miscut angles /? from —0.30° to 0.30°. It is 
noteworthy that a tiny variation of miscut angles ft has no effect 
on the values of K 2 and iC u (Figure 6(a)), whereas significantly affects 
the fitted values of the first-order magnetocrystalline anisotropy 
^(Figure 6(b)). Since the global AMR properties are measured for 
the whole sample, the use of the mean miscut angle of Si(lll) of 
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Figure 5 | The normalized magnetic torque curves. Normalized magnetic 
torque curves superposed by the first- and second-order 
magnetocrystalline anisotropics and uniaxial anisotropy at different fields 
(a) and the normalized torque contributed only by the first- and second- 
order magnetocrystalline anisotropics at field of 1600 A/m (b), the solid 
lines are fitting curves. For comparison, the normalized magnetic torque 
curves were also plotted in Figure 5(a). 

about 0.10° is reasonable. The fitted value of K Y = 3.4 X 10 4 J/m 3 is 
comparable with the value of K Y = 4.5 X 10 4 J/m 3 for bulk bcc-Ve. 

The AMR results are consistent with the fitting results from FMR, 
where K x = 4.4 X 10 4 J/m 3 , K 2 = 2.2 X 10 3 J/m 3 and K u = -5.9 X 
10 2 J/m 3 for 45 ML Fe film on Si(lll) 16 . For comparison, we also 
measured the 45 ML Fe sample in Ref. 16 by FMR to cross-check the 
measured anisotropy constant obtained by AMR. As shown in 
figure 5(a), the magnetic anisotropy constants K Y = 4.4 X 10 4 J/ 
m 3 ,K 2 = 2.1 X 10 3 J/m 3 and K u = -1.1 X 10 3 J/m 3 can be obtained. 
Both AMR and FMR give almost same values of K x and K 2y while the 
absolute value of K u obtained by AMR are slightly larger. A careful 
comparison between these two techniques to measure the magnetic 
anisotropy constants is in progress. 

The negative value of K u suggests that the easy axis is perpendic- 
ular to the step direction, which is in good agreement with the hys- 
teresis loops (Fig. 3 (a)). By using the magnetic anisotropy constants, 
the angular dependence of coercivity and remenance were also simu- 
lated in term of coherent rotation magnetization reversal, as illu- 
strated in Figs. 3(b) and (c), respectively. We can observe that the 
calculated remanence is consistent with experimental one, whereas 
the calculated coercivity deviates the experimental one significantly 
in the easy axis direction. The significant deviation implies that the 
magnetization reversal is not governed by coherent rotation, and 
consequently, magnetic hysteresis loop cannot provide a detailed 
symmetry of magnetic anisotropy energies. 

The contribution of K Y and K 2 to magnetic torque is about one 
order of magnitude smaller than that of uniaxial magnetic anisotropy 
constant K u . In order to separate their contributions, Fourier analysis 
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Figure 6 | The fitted magnetic anisotropy constants K 2 and K u (a), and K x 
(b) for various miscut angles /?from —0.30° to 0.30°. 

and inverse Fourier transform are used to analyze the torque curve. If 
we deduct the contribution of the uniaxial magnetic anisotropy con- 
stant K u to the magnetic torque, it is obviously shown in figure 5(b) 
that the normalized magnetic torque curve Ik u k 2 (!?m) * s tne super- 
position of a four-fold and a six-fold anisotropy contributed only 
from first- and second-order magnetocrystalline anisotropy con- 
stants, respectively. 

Discussion 

We present the MOKE and AMR measurements of Fe single crystal 
film on 0.1° vicinal Si(lll) substrate with atomically flat ultrathin 
p(2 X 2) iron silicide as buffer layer. Unfortunately, owing to the 
small coercivity (< 10 Oe) in the sample, only the two-fold symmetry 
in remanence and coercivity can be observed, while the contribution 
of magnetocrystalline anisotropy constants K Y and K 2 cannot be 
determined from MOKE measurements. On the other hand, the 
AMR results demonstrate that the symmetry of magnetocrystalline 
anisotropy energy changes from the six-fold to a superposition of six- 
fold, four-fold and a weakly uniaxial contribution due to the tiny 
misorientation from Fe(l 1 1) plane. Although the use of AMR to mea- 
sure the magnetic anisotropy was introduced long time ago 20 " 26 , to 
our knowledge, a precise determination of various magnetic aniso- 
tropy constants of Fe(l 1 1) film on Si(l 1 1)-7 X 7 surface with so small 
miscut angle was not reported in literature. The fitted value of the 
first-order magnetocrystalline anisotropy K Y is significantly influ- 
enced by the tiny variation of miscut angles On the other hand, 
the values of K 2 and K u are unchanged. Our work suggests that the 
AMR measurements can precisely separate the detailed contributions 
of various magnetic anisotropy constants of single crystalline Fe ultra- 
thin film grown on vicinal Si(lll) surface. 

Methods 

The sample was prepared on Si( 1 1 1 ) wafers with nominally miscut angle of 0. 1 0 along 
[11-2] using an ultrahigh vacuum molecular beam epitaxial chamber (MBE) 
equipped with the scanning tunneling microscope(STM) and low-energy electron 
diffraction(LEED). The base pressure of MBE is kept around 2 X 10" 10 mbar and all 
the experiments were conducted at room temperature. After a well-established pro- 
cedure 1516 , the well-defined reconstructed Si(lll)-7 X 7 surface was obtained. The 
buffer layer was deposited on the wafers for 1.5 ML of Fe (99.999% purity) heated by 
e-beam bombardment with a deposition rate of 1.5 ML/min, then annealed at 700 K 
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for 10 min. This procedure gives a highly ordered 2X2 periodic iron silicide 
structure to prevent the Fe/Si intermixing 6 . Fe film with thickness of 21 ML was 
deposited on the iron silicide template. The STM (VT-STM) measurements were 
performed at the Si substrate, iron silicide template and the Fe film. A non-magnetic 
NaCl with thickness of 14 ML was deposited on the sample as a capping layer to 
protect samples oxidization. 

The MOKE measurements were carried out at room temperature and described in 
detail elsewhere 1516 . The home-made AMR setup consists of a Wheatstone bridge, a 
Lock- in amplifier (Stanford Research Systems SR830 DSP), a temperature controller 
(Stability < 0.0012°C/h), and a rotational sample stage. Magnetic field is provided by 
a Helmholtz coil. In the experiments, a sufficiently large and stable field is applied to 
guarantee a true single-domain behavior of the specimen. The application of 
Wheatstone bridge and highly stable temperature controller ensures the sensitivity of 
AMR better than 0.01% in the entire measurements. The sample size is 3 mm X 
5 mm for AMR measurements, which were performed with a standard four-point 
method. 
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